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1 Introduction 

For higher dimensional considerations of gravity, Lovelock generalization is the most natural 
and appropriate as equation of motion continues to remain second order despite Lagrangian 
being homogeneous polynomial in Riemann. This property singles out Lovelock from all 
other generalizations and it includes Einstein gravity in linear order = 1 where N is 
order of Lovelock polynomial action. It is well known that Einstein gravity is kinematic in 
d = 2N + 1 = 3 dimension which means Ricci equal to zero implies Riemann zero. In other 
words, Riemann is entirely determined in terms of Ricci. Could this property be true for 
all N]i.e. does there exist a higher order generalization of Riemann such that in all odd 
d = 2N + 1 dimensions it is entirely determined by its Ricci? The answer is yes [1] for 
pure Lovelock gravity for which the action contains only one A^th order term, besides of 
course the cosmological constant A, in the action (there is no sum over N). We define [2] 
A^th order analogue of Riemann (trace of its Bianchi derivative yields analogue of Einstein 
tensor and agrees with the one that is obtained by varying A^th order Lovelock action) which 
is a polynomial in Riemann and its contraction giving corresponding Ricci, and so we have 
Lovelock Riemann and Ricci tensors. Then it turns out that whenever Lovelock Ricci vanishes 
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in odd d = 2N + 1 dimensions, so does Lovelock Riemann. That is, pure Lovelock gravity 
which includes Einstein gravity for = 1, is kinematic in all odd d = 2A^ + 1 dimensions. It is 
Lovelock kinematic because it is Lovelock Riemann that vanishes and not the usual Riemann. 
Kinematicity so defined is thus a universal gravitational property in higher dimensions, and 
if it is to be respected always, for a given N there are only two pertinent dimensions d = 
2N + 1,2A'" + 2. For d < 2N + 1, it is vacuous while for d > 2N + 2, it is next order 
(N + l)th that should be considered else kinematicity property will be violated. We shall 
therefore adhere to universality of this property and shall emphasize the pure Lovelock case 
for d = 2A^ + 1, 2A^ + 2. 

Spherically symmetric static black holes in Gauss-Bonnet (GB) gravity have been studied 
[3, 4] and in fact for Lovelock gravity in general [5, 6]. It turns out that ultimately one needs 
only to solve an algebraic polynomial of degree N. To manage this polynomial and for 
other physical considerations, it was assumed [7, 8] that all coupling constants are given in 
terms of the unique ground state, A. This made the equation degenerate and hence could 
be trivially solved. On the other hand, if we consider pure Lovelock equation which has only 
one A^th order term with A, then the equation is derivative degenerate and again could be 
solved trivially [9]. In all this, black hole horizon was a sphere of constant curvature. Some 
authors have considered more general cases with arbitrary couplings for each Lovelock term, 
see [10, 11], though an exact solution for polynomial equation cannot be obtained beyond 
A^ = 4. 

The next generalization came through consideration of horizon topology being a prod¬ 
uct of two spheres. Now horizon space is non-maximally symmetric Einstein space which 
has non-zero Weyl curvature. However both Weyl and Riemann have vanishing covariant 
derivative. Its first example was given by Nariai metric [12]. An interesting GB black hole 
solution with two sphere topology was obtained by Dotti and Gleiser [13] and subsequently 
its various features like uniqueness and stability were studied by several authors [14-17]. It 
turns out that non-zero Weyl curvature makes a non-trivial contribution in solution in terms 
of a constant which is negative for GB case. This implies that A should always be present and 
positive and there also occurs non central singularity. To avoid that, a range is prescribed for 
black hole mass in terms of A [18]. 

In this paper our main aim is to extend this framework to pure Lovelock of arbitrary 
order. As a matter of fact, our study goes beyond pure Lovelock case {d = 2N + 2) with 
horizon topology x to include any d = 2do + 2, do > N with horizon topology 
S(do) 

X That is, we solve Lovelock A-vacuum equation of arbitrary order N for a static 

black hole with two sphere topology in a spacetime with even dimension d > 2N + 2. It is 
remarkable that algebraic character of ultimate equation is carried through for this enlarged 
framework as well. The product of two spheres produces solid angle deficits which cancel out 
each-other’s contribution in Ricci tensor but not in Riemann. For A^ > 1, there is always 
Riemann present in the equation and to cancel out solid angle deficit contribution to it, a 
constant, say p, is required which is zero for Einstein N = 1 case. Further it turns out that 
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A and p always bear opposite sign, the former is negative and positive according to N being 
even and odd, and vice-versa for p. Whenever N is even, p < 0, A must be present and 
positive. To avoid non central singularity and for existence of horizon, black hole mass has 
to lie in a range prescribed in terms of A. On the other hand for odd N when p > 0, there 
is no such constraint. It is remarkable that the similar discerning feature also emerges from 
thermodynamical stability considerations. It favors odd N and A < 0 as against even N and 
A > 0. It is remarkable that two sphere topology distinguishes between odd and even N and 
consequently negative and positive A quite clearly. This is a new realization that has come 
up only when we have considered both odd and even N > 1. 

The paper is organized as follows: In next Section we set up equation of motion for 
pure Lovelock A-vacuum of arbitrary order N for a static metric with two sphere topology, 
and obtain general black hole solution. It is followed by consideration of particular cases for 
N = 1,2,3,4. In Sec. IV we compare two sphere black holes with the corresponding one 
sphere ones while Sec. V is devoted to thermodynamical considerations. We end up with 
discussion. 

2 Equation of motion 

Consider static metric configuration with product horizon space, x having sym¬ 

metry group SO (do -|- 1) X SO(do 1), 

ds^ = -A{r) dt^ + ^ dr^ + , (2.1) 

describing an even-dimensional spacetime with dimension d = 2do-\-2, where is metric of 

do dimensional unit sphere. We compute Riemann tensor for the metric and its nonvanishing 
components are as follows: 

Di0,l) = = -\A"{r) 

D{0,i) = Ror = -^A'ir) 

2 r 

D{0,i') = Ro^^^' = -^A\r) 

2r 

Dil,i) = Ril^ = ~A'{r) 

D{l,i') = Rii^^' = -^A'{r) 

2r 

D{i,j) = mA = ^{1-A{r)) 

= = ^{1-Air)) 

Dii,i') = R^r' =-^Air) (2.2) 
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where indices 0,1 refer to coordinates t, r, and indices i, i' to angle coordinates of two spheres. 
Notice that there is a crucial difference, in fact the only one, between the one-sphere case, 
dS‘^d- 2 )^ and our present two-spheres case, which is the last entry, D{i, i'), taking indices from 
two spheres. All the other components are the same for both configurations. 

For this pure Lovelock Lagrangian of order N, we will show that the equation of motion 
(EOM) for the metric (2.1) is ultimately given by a single algebraic equation. Lovelock 
Lagrangian of arbitrary order N is given by 

= , (2.3) 

which leads to EOM 

[R-Nrp = 5;) A = o, (2.4) 

where A is, up to some combinatorial factor, the cosmological constant. We keep the factor 
for reasons of later convenience. 

Our aim is to implement Eq. (2.2) into EOM (2.4) to get a differential equation for A{r). 
It is obvious that for non-vacuous EOM, fJ- = p, and we thus write =: [Rn]p,- From Eq. 

(2.2), it is clear that 

= [^iv]o (2.5) 

whereas for angles p = i^i' 

[CN]i = [CN]i' =: [^iv]i ,Vi, i'. (2.6) 

So we are left with two equations, [£Ar]o = 0 and [C-n]^ = 0. With some manipulations, which 
we leave for Appendices A and B, it can be shown that the former is in fact a total derivative 
(Appendix A) while the latter takes the form (Appendix B), 

[Rn], = ^dr(^[CN]oy (2.7) 

An equation of this kind is indeed expected on the basis of Noether identities for gauge 
theories. Note that here the gauge freedom is general covariance. In the case N = 1 Noether 
identity Eq.(2.7) is the contracted Bianchi identity for Einstein tensor. Thus we need to 
consider only the equation [TatJo = 0 . 

What is truly noteworthy is that [TtvJo is a total derivative. It is well known that this is 
already the case with one sphere topology [5, 6]. It is remarkable that this interesting feature 
carries forward also to the present case of two sphere topology (See Appendix A). Now the 
A-free term in Eq.(2.4) for [TatJo has the form ^ ^r^'^°“^'^+^PAr(A(r))^, where PN{A{r)) is 
a polynomial of degree N. Here the factor has its origin in the density y/—g. Then the 
equation, [TatJo = 0 from Eq. (2.4) becomes 
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which trivially integrates (with the appropriate redefinition A = — -) to give 

2oo + 1 

P«(^(r)) = Ar='' + ;35^, (2,9) 

where M is an integration constant. Note that in asymptotic expansion, the potential due 
to M will go as which is the Schwarzschild potential in d dimension. 

This formally indicates that M is mass of the configuration. This is a general feature of pure 
Lovelock gravity that even though EOM is free of Einstein term yet the potential due to mass 
approaches asymptotically the Schwarzschild potential in the corresponding dimension [19]. 

As emphasized earlier, it is very remarkable that the general solution can now be obtained 
by simply solving the above algebraic equation of degree N. 


The polynomial Pi\f{A{r)) takes the following form, 

PNiA) := ^ C[N - do - 1] (-^)'(1 - 


N-l 


1=0 


( 2 . 10 ) 


with the combinatorial object 

- m 

C[m, s] = —Y, 


m 


s\ 


Si 


2"^ \ k J {s — 2m + 2 k)l {s — 2k)\ 


( 2 . 11 ) 


(An extension of the factorial to negative numbers is understood in this formula (2.11), using 
the gamma function. In particular notice that C[Ai, Al] vanishes for odd N) 

This polynomial can also be expressed alternatively as 


N N-l 


Pn{A) = 

1=0 k=0 


1 A^! 

2^ k\l\{-k-l + N)\ 


_ 

(do “ 2/c 


do! 

0! {dQ + 2k-2{N -1) -l)\ 




( 2 . 12 ) 


Let us write 


and define Qn{^) 


A{r) 


dp — 1 
2do - 1 


(l-vl/(r)), 


(2.13) 


Pj\f{A). With this substitution the algebraic equation (2.9) becomes 


A/f 

Qn{^) := , (2.14) 

The substitution (2.13) has the specific property that the polynomial (5Ar(T) is now free of 
term. This means that sum of its N roots T;, / = 1, • • • Ai is zero. 

In terms of T, the metric is written as 


ds^ = -(1 - 'l'(r)) dt^ + ■ (2-15) 
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Eq.(2.14) is the general algebraic equation for pure Lovelock black hole of arbitrary order 
N. Next we consider specific cases for N = 1,2, 3,4 to show that the known Einstein and 
Gauss-Bonnet black holes [3, 4] are all included in this equation. 


3 Particular cases 

3.1 Case N = 1: Einstein black hole 

In this case, Eq. (2.8) takes the form 


giving 

, , , 9 M , , 

( 3 - 1 ) 

for the metric (2.15). This is the usual Schwarzschild-dS/AdS solution except for constant 
factor 2d(~-i before two spheres metric. Note that this factor causes solid angle dehcit for 
each sphere but the two together conspire to cancel out each-other in Ricci giving rise to 
A-vacuum. This feature was already noticed by [20]. However this cancellation cannot carry 
forward to higher order Lovelock simply because analogue of Ricci or Einstein tensor for 
higher order Lovelock would also involve Riemann. Their contributions to Riemann do not 
however cancel out each-other. For that an additional constant would be required for > 1 
solutions as as we would see in the cases that follow. 


3.2 Case N = 2: Gauss-Bonnet black hole 


In Gauss-Bonnet case, the substitution (2.13) turns out to be very efficient, making the 

' Then the solution for Eq. 


polynomial (2.10) proportional to ^'^(r) -|- 
(2.8) is of the form 


(do-1)2(24-3)' 


'l'(r) = Ayp -|- Ar^ -|- 


M 


j^2dQ—S 


where 


P = 


do 


(3.2) 


(3.3) 


{do - 1)2(24 - 3) 

We choose the upper positive sign in solution (3.2) for existence of black hole horizon as well 
as for gravity being attractive. Note here the constant p under the radical is negative and it 
is the one that compensates the term coming from square of Riemann in GB Ricci analogue 
which is quadratic in Ricci and Riemann. As argued earlier for Einstein case that solid angle 
deficits of the two spheres cancel out each other to leave A-vacuum intact while in the case 
of GB, there is also a quadratic term in Riemann, it is that which requires this additional 
constant. For pure GB black hole, we have d = 2N -|- 2 = 6, 4 = 2 and p = —2. This case 
has been comprehensively discussed in [18]. That is A and p would always bear opposite sign 
and they are respectively positive and negative for even and odd N and vice-versa. We shall 
verify this feature in the following cases for = 3,4. 


- 6 - 








3.3 Case N = 3: Third order pure Lovelock black hole 

With the appropriate redefinitions, always preserving the signs of A and M, Eq. (2.9) is now, 


with 


a + 6 'l'(r) — p = Ar® + 


M 

j^‘2idQ —5 


(3.4) 


a= {do- lf{2do - 5){2do - 3) > 0 , 
b = 3{do - l)do(2do -5) >0, 

p = 2do(2do + l) >0. (3.5) 


The critical dimension (pure Lovelock) for N = 3 is d = 2N + 2 = 8, do = 3, and so we write 
(after redefining A and M to absorb a factor of 6) 


4T(r)^ + 3T(r) 


7 + Ar^ + 


M 

r 


For existence of horizon, '^{vh) = 1, we need Ar® + ^ = 0 which gives M = —A.rJ^. The 
equation then becomes 

4'I'(r)^ + 3'I'(r) = 7+— (r’^ — rl). (3.6) 

r 

It is clear that + 3 4' is monotonically increasing. To keep T < 1 (i.e, A{r) > 0), which 

Ay 7 

in Eq.(3.6) means RHS < 7, we need —{r — rL) < 0. For r > rh, this means that A must be 

r 

{ M \\ 

negative A < 0 which also ensures M > 0 in = I —- I . The solution T (r) takes the form 


'l'(r) 


(yWri + Sy^"^ -1 


with 5 = 7 — ■^(r'^ — r^) . Note that r > r/i ^ 5 < 7 4'(r) < 1. This describes the region 

outside black hole horizon. Thus there is no constraint on horizon radius and mass is given 
by M = |A|r^. 

The limit r/j —)• 0 (=► M —)■ 0) is regular for 'l'(r) but then the background develops a 
naked singularity at r = 0. In fact the Kretschmann scalar, K = , behaves for 

r —>■ 0 as 

24 

(3.7) 

Note also that in the r —>■ oo limit, 4'(r) —>■ — '^^ 2/3 mentioned earlier, p = 7 > 0 and 

A is required to be negative. On the other hand for N = 2 GB case, it was the reverse, 
p = —2 < 0 and A > 0. That is, A has always to be non-zero, positive for even N and 
negative for odd N while for = 1 it could be both positive and negative. Pure Lovelock 
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critical dimension for = 1 is d = 4 for which do = 1 that makes x 5(1) space flat. 
Therefore dimension has always to be > 4 and hence Einstein case in this setting is not a 
proper pure Lovelock critical dimension case. We will also see in the next case of At = 4 that 
A > 0 and p < 0. 


3.3.1 The noncritical case do > 3 

The N = 3 case, for do > 3 (non pure Lovelock) allows for a window of positive A. Although 
it is not the case that really interests us (which is do = At = 3), we will discuss it here for the 
shake of completeness. 

Let us assume do = 4, which is representative of the general do > 3 case. Equation (3.4) 
is (an irrelevant positive factor has been eliminated from the dehnition (2.14)) 

Q,{^{r))=Ar^ + ^ (3.8) 

What happened in the critical, pure Lovelock case was that at the horizon ('I'(r/i) = 1) the 
polynomial Qs vanished (the reason is that Pn{0) = 0 for odd do = At). Now instead we have 
< 53 ( 1 ) > 0 . and the polynomial is monotonically increasing, Q 3 ('I') > 0 . 

At the horizon we have Q 3 (l) (= ^* 3 ( 0 )) = Ar® + ^. Thus we can get M = P 3 ( 0 )r)) — Ar® 
and so we write (3.8)) as 

Q 3 ('k(r)) = A(r® - i) + ^ 3 ( 0 )^ . (3.9) 


Since 'I'(r/i) = 1, for r > but close enough to the horizon, we need 4'(r) < 1, otherwise 
A(r) turns negative. Using that the polynomial (53('L) is monotonically increasing, this means 
that we need, at the horizon, J;Q 3 (^(’’))|r^ < 0 which, in terms of the rhs of (3.9), yields the 
condition 


A < 


^ 3 ( 0 ) 

34 


Thus there is a window for positive A. Of course any negative value for A is admissible and 
in such case the analysis continues along the same path as the critical do = N = 3 case. The 
case with positive A can be parametrized by taking A = with o > 0. We then find 

that there is a minimum rg for 'l'(r) ('I''(ro) = 0) at rg = (1 + |)9r/i > r^. Then, for r > rg 
'I''(r) turns positive and we obtain a second, cosmological deSitter-like horizon at Vc when 
4'(rc) = 1. 


3.4 Case N = 4: Fourth order pure Lovelock black hole 

Again with the appropriate redefinitions of A and M, Eq. (2.9) takes the form. 


a'I'(r)^ + 5T(r)^ + c'I'(r) — p = Ar® + 
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Figure 1. Plot of F(4') (Eq. (3.11)). 

Horizontal lines are y = F'(l) = 37 and y = F'('I'o) ~ —7.81 at the minimum 'I'o — 0.41. 


with 


a = (4-1)^(2do-7)(2do-5)(2do-3) >0, 
b = 6(do - l)^do(2do - 7)(2 do - 5) > 0 , 
c = -8(do-l)do(2do-7)(2do + l) <0, 

p = -3do(do(10do + 3) + 2) <0. (3.10) 


All cases for do > 4 are similar and we focus on the pure Lovelock N = 4, do = 4 case. 
Redefining again for convenience A and M, the equation for 'l'(r) becomes 

239 A M 

F(4'(r)) = 45^(r)^ + 24^(r)2 - 32^(r) =-+ -r® + — . (3.11) 

3 8 r 


The plot of F{^) against T (see Fig. 1.) has a parabola like shape with a minimum at 
'I'o — 0.41, where F{^o) takes the approximate value F{^o) ~ —7.81. We need to find a 
region for r for which 4'(r) < 1. Note that T(l) = 37. Thus we need values of r such that, 
approximately 


or 


232 A a M 

-7.81 <-+ + — < 37, 

3 8 r 


A c M , , 

69.2 < -r®+ —< 114.3. (3.12) 

Clearly this inequality requires A > 0 while M is always to be positive. Otherwise one of 
the bounds on r will be an unwanted naked singularity, 'l'(r) turning complex. In fact the 

minimum for r should respect the range given in (3.12). The minimum is at tq = ^x) 
requiring that it satisfies the bounds (3.12), we end up with following range for A and M, 


61.8 < mIas < 101 ,6. 


(3.13) 


- 9 - 







For A and M satisfying this bound, we can find a region < r < rc such that 'h(r/j) = 
'h(rc) = 1 where is black hole and Tc cosmological horizon (similar to deSitter), and tq also 
satisfies < vq < r^. Note that here p = —232/3 and A > 0 bearing out the general pattern 
as envisaged earlier. 

3.5 The general pattern 

As argued in Introduction, the distinguishing property of pure Lovelock gravity is its kine- 
maticity in odd d + 2N + 1 dimension. For this property to be universal;i.e. true for all odd 
d = 2N + 1 dimensions, the equation should be valid only in the two critical odd and even 
dimensions, d = 2N + 1,2A^ + 2. We shall therefore stick to the critical, pure Lovelock case, 
with do = N. Based on the cases Ai = 2,3,4 and more (we have explored up to = 12 but 
nothing prevents from going to any higher N), there emerges the following general pattern 
that A and p (note that up to a positive factor which will be absorbed in A and M, p is 
defined as p = —Qn{0)) always bear opposite sign, the former is respectively positive and 
negative for even and odd values of N while the opposite is true for the latter. For N even 
the shape of ( 57 v('k) is parabola with a single minimum 'km satisfying 0 < 'km < 1, similar 
to the case A^ = 4 as shown in in Fig. 1. For N odd Qn{^) is a monotonically increasing 
function. 

It may be noted that the parameter p arises from countering the deficit angle contribution 
in the Riemann tensor and hence its sign and value is determined by this geometric feature. 
On the other hand for even Al, A > 0 is required for the polynomial (2.14) to have a minimum 
while for odd N, for some r, (5Ar('k) < C[N, A^] = 0 in the critical do = N dimension, and for 
which A has to be negative. 

As expected whenever A > 0 {N even), there exist both black hole and cosmological 
deSitter-like horizons, and a range of values for black hole mass in terms of A. When A < 0 
{N odd), there is only black hole horizon without any constraint on mass (See Sec. 5 for 
non pure Lovelock cases). As a matter of fact, since C[N,N] = 0 for N odd, we will have 
Q7 v(1 ) = 0 for odd do = N because (57v(l) = Pn{0) = C[N,do]- Thus it is both BH 
and cosmological horizons in d = 6,10,14,... and only one BH horizon in d = 8,12,16,... 
dimensions. 

4 Comparison with one sphere case 

It is well known that static spherically symmetric vacuum equation for the general Lovelock 
Lagrangian ultimately reduces to an algebraic equation [5, 6]. The difference between one 
and two spheres configurations is the absence of the terms with i! indices for the former 
in (2.2). This fact greatly simplifies the analysis and it follows on the same lines as in 
our case. The analogue of of the polynomial (2.10) becomes just a monomial in (1 — A), 
namely P/v(A)T) ;= — A)^. The combinatorial factor is customarily absorbed 

in the general Lovelock Lagrangian £lov. = with a redefinition of coefficients cn- 
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(4.1) 


Writing A{r) = 1 — r‘^B{r), we have the analogue of EOM (2.8) as 

^|'^2do+l^(^)iV^ ^^2doA. 


No 

If the Lagrangian is a sum of several Lovelock terms E catEat, then B{r)^ in (4.1) is 

N=1 

substituted by YlcNB{r)^ and we have, on integrating (4.1), 

E (4.2) 

N=0 


where the coefficient cq is proportional to —A and the integration constant M is proportional 
to black hole mass. 

This nice simplification in Eq. (4.2), due to the fact that the power of the variable r in 
the bracket in the Ihs of (4.1) does not depend on N, is no longer available in the two-spheres 
case because of the different structure of the polynomial. The crucial difference in the two 
cases is the presence of terms of the type D{i,i') in Eq.(2.2) which, unlike the terms D(i,i) 
and are not proportional to 1 — A(r), but to A{r). 


5 Thermodynamics 

Even without knowing solution to Eq. (2.9) explicitly, it is possible nevertheless to perform 
thermodynamical analysis for its black hole (BH) solutions by just using EOM relying on the 
fact that these equations are of polynomial form. Therefore the method described below can 
be applied whenever EOM are of this type. Here we assume that solution A{r) exist with a 
horizon A{rh) = 0, for some range of parameters A and M. As a matter of fact we already 
know that such BHs exist for = 1, 2 and we have shown in Sec. 3 that this is indeed the 
case for arbitrary N. 

Let us consider Eq. (2.9) for a putative solution A{r, M) (since A plays no special role 
here we do not need to explicitate the dependence of A on A). Thus we have the identity 

= + (5.1) 

and next we differentiate this identity w.r.t. M and write We get 

, dA 1 

^nW = ^2do-2N+l ■ 

We also note that ^^(0) < 0 because the definition (2.10) gives 

P^(0) = -N {dl C[N - 1, do - 1] + C[iV, do]) < 0. 


dA 


A=0 — — 


1 


Thus 




(5.2) 








On the other hand, we can trade and M for each other, which defines the function M{rfi) 
so that A{r,M{r)) = 0 while varying r. Thus we can write the identity 


A{r,M{r)) = 0 


as implicitly defining the function M{r). Its r-derivative gives the identity 


BA 

A'(r,M(r)) + —U=oM'(r)=0, 


where T'(r, M{r)) is the derivative with respect to the first argument r. In fact the Euclidean 
method identifies A'(r,M{r)) as the Hawking temperature A'(r,M{r)) = 47rT(r). Thus we 
have 


M'(r) 


= —dvr 


1 


dA 

dM 


A=0 




(5.3) 


5.1 Entropy 

The entropy is then obtained by integration of the First Law 

fdM^ T'* ^ ^ 47r|P;,(0)| odo-27V+2 

J T Jo T{r) 2do-2N + 2 ^ 


(5.4) 


where we have assumed the vanishing of the entropy when the horizon shrinks to zero. 

Since we only know that M is proportional to the mass of the BH, there is an undeter¬ 
mined numerical factor not included in (5.4). The relevant finding nonetheless is that S is 
proportional to ^‘^o-‘iN+ 2 ^ Note that in the critical Lovelock dimension, d = 2N + 2, the 
entropy becomes proportional to or, since the area A of the horizon is proportional to 
, it turns out that 


S, 


(do^N) 


Av , 


which confirms that we are in the same universality class as in the spherical pure Lovelock 
BH, as discussed in [9]. 


5.2 Temperature 

The temperature con be obtained using the identity (5.1), in which M is a fixed parameter. 
The r-derivative of (5.1), written at the horizon rh gives 

P^(0) A'irf,, M) = 2A Arf-1 - {2do -2N + 1) 2 do^7V+2 > 

'’'h 

which means that, since Pj>^{0) < 0 and A'[r, M{r)) = 47rT(r), 

^ ^ 47r|P;(,(0)| + l)^Mo-2iV+2 - 2AAr2^ . 
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We can eliminate A from this expression by taking into account that (5.1) implies, at the 
horizon, 

Pn{0) = + 2do-2N+l ’ 

(Note incidentally that M{rh) is trivially obtained from (5.5). Note also that the value 
PNiO) = C[N,do] is obtained from (2.10).) 

Using (5.5) we obtain 


4vr|PA,(0)| 


M 


2do-2N+2 



(5.6) 


which in the critical Lovelock dimension, d = 2N + 2, gives 

1 


X 


(dQ=N) 


4vr|P^(0)| 


2N ^ 

2N + 1)—^ — 


h 


Th 


which agrees, up to numerical factors, with the results obtained in [9] for the spherical pure 
Lovelock BH. 


5.3 Stability 

Here we follow the approach of Refs [11, 17]. Local thermodynamical stability will correspond 
■n . ^ dM , . 

to positive specihc heat, Cg = compute it as 

_ 

" T'{rn) ■ 

Using Eq. (5.5) to express M in terms of r^, we have 

Mir^) = {Pn{0) - (5.7) 

Then temperature in Eq. (5.6) is given by 

To keep M > 0, as required from the analysis in Sec. 3, 

P^(0)>Ar2^. (5.9) 


Note that PAr(O) > 0 always. It only vanishes at the critical do = N dimension for odd N. In 
all other cases Pn{0) > 0. To keep T > 0, we must have 


Pv(0) > 


(2do + 1) 

{2do -2N + 1) 




.2N 


(5.10) 
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Thus for A < 0, the above constraints (5.9) and (5.10) impose no restrictions whereas for 
A > 0 the constraint (5.10), which anyway implies (5.9), must be satisfied. 

From Eq. (5.3) it is clear that T > 0 M' > 0. To have local thermodynamical stability 
we therefore need T'{rh) > 0. In fact, since 

( - <2* - ^ - <2* + 1)(27V - 1)A rf-2) , 

the condition T'{rfi) > 0 is equivalent to 


Pn{0) < 


{2do + 1)(2A^ — 1) 2N 

2do-2N + l ^ 


(5.11) 


Local thermodynamical stability is thus guaranteed if the bounds (5.10) and (5.11) are sat¬ 
isfied. Since Pn{0) > 0, Eq.(5.11) is inconsistent for positive A > 0, and hence it is ruled 
out. Thus we need A < 0;i.e. N odd. If in addition do = N, PAr(O) = 0 and therefore no 
restrictions of any kind. If FW(0) > 0, there is a lower bound (Eq.(5.11)) on BH size. 


rr> 


(2do + l)(2A^-l)|A| 


Ejv(O) 


(5.12) 


consequently its mass, as given in Eq.(5.7), is also bounded from below. 


Next we consider global thermodynamical stability by requiring free energy to be negative. 
Using Eqs (5.4) and (5.8), we compute free energy, 

F = M-TS= ^ ^2Ar - l)Arfo+i) . (5.13) 

Eor do = N, Pn{0) = 0 and A < 0 {N odd), hence free energy is always negative ensuring 
global stability while for even N and A > 0 it is always positive thereby implying instability. 
That is black hole is thermodynamically unstable for A > 0 and N even. Thus thermody¬ 
namical stability requires A < 0 and N odd. When do = N (pure Lovelock), there is no 
restriction of any kind while for do > N for which F’Ar(O) > 0, there is lower bound on horizon 
size given by 


> 


Pn{0) 


(5.14) 


(2A^- 1)|A| ’ 

which already implies, as expected, local thermodynamical stability Eq.(5.12). 

Note however that our analysis is formal in the following sense. When do > N the 
configuration with M = 0 is not allowed in our formalism as it is forbidden from equations 
(5.7) and (5.14), since F7v(0) > 0 . One could think that in the critical, odd do = N case, the 
limit —7- 0 which sends also M to zero, could be taken without restrictions. And in fact 
the entropy (5.4) and the temperature (5.8) (with Pn{^) = 0) vanish in this limit. But then 
a naked singularity appears at the origin r = 0 (see eq.(3.7) in section 3.3) and thus we must 
keep M > 0 to prevent this singularity from forming. This means that the global stability 
analysis must be considered formal. 
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6 Discussion 


We have shown, within the framework of Lovelock theories, that the property that the EOM 
are polynomial, which it is well known to hold in the case of spherically symmetric horizons, 
can be generalized to horizons that are product of two equal spheres (though now the Weyl 
tensor is not vanishing), which are a particular case of Einstein manifolds. On the other 
hand we also know that in the case of Gauss-Bonnet gravity (Lovelock N=2) this polynomial 
property of the EOM holds whenever the horizon is an Einstein manifold [13] (thus including 
the one-sphere, two-spheres cases). It can be conjectured that, in a general Lovelock theory, 
the polynomial property that we find for two-spheres topology could be extended to the more 
general case of Einstein horizons. An indication in this direction is the study of = 3 
Lovelock in [21]. 

What does two sphere topology entail is that horizon space is non-maximally symmetric 
Einstein space which has non-zero Weyl curvature, though both Weyl and Riemann have 
vanishing covariant derivative. Product of spheres cause solid angle deficits which cancel out 
each-other’s contribution to Ricci tensor [20] but not to Riemann tensor. Thus for N > 2, 
EOM always has Riemann present and hence contribution to it from solid angle deficit has 
to be anulled out, and that is done with an additional numerical parameter p, as shown 
in the N = 2 case, in which it appears under the radical in the solution. It is in fact a 
measure of non-zero Weyl curvature [13] which is absent for = 1. The general feature 
of vacuum equation, ultimately reducing to an algebraic polynomial is carried through even 
for two sphere topology for do > N. That is why we have been able to solve pure Lovelock 
vacuum equation for arbitrary order N in this enlarged general setting. Though our general 
analysis refers to the critical dimension do = N, however the equation reducing to an algebraic 
polynomial is a general feature true for all do > N. 

It turns out that in the critical do = N case, constant p and A always bear opposite sign, 
the former is positive or negative according as N is odd or even and opposite is the case for A. 
For even N, p < 0 and A > 0, spacetime has two horizons, black hole and cosmological, and 
black hole mass is constrained in a range given in terms of A. While for odd Al, A < 0, there is 
only one black hole horizon without any restriction. Remarkably thermodynamical stability 
also has the similar discerning feature. That is spacetime is thermodynamically stable for 
odd N and A < 0 while it is unstable for even N and A > 0. 

We thus conclude that two sphere topology for pure Lovelock gravity clearly discerns 
between odd and even Lovelock order N, consequently negative and positive A and it is the 
former which is favored. This is a remarkable and interesting property of this setting which 
has been revealed only when we probed higher orders, N > 2. Let’s reiterate that universality 
of kinematic property of gravity in all odd d = 2N +1 dimensions demands the pure Lovelock 
equation which is valid only for the two d = 2N + 1, 2A1-|- 2 dimensions. This is the main and 
critical motivation for pure Lovelock gravity. In the context of the present study we therefore 
always take d = 2A^ -|- 2 which implies do = N. 
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Appendix A: Proof that [£Ar]o is a derivative 


Consider equation (2.4) for p = /i = 0. The factor is, up to irrelevant dependences on 
the angles, . Obviously the term with the cosmological constant is an r-derivative, so we 
will only consider in this appendix the expression [TatJo with vanishing A, 

[-C-ATJO — ' O0pi...pjvcri...o-jv ■ ■ ■ ^flNVN ■ 

This expresion is a sum of products of N components of the Riemann tensor. We must fill in 
it the components (2.2). Since the index 0 will not be present in these components, we can 
only introduce the factors Note that the first two 

objects in this list contain A'{r), whereas the rest only contain A{r). So (6.1) will become 
the sum of two objects, one of them containing A'{r) as a factor and the other without this 
derivative. To construct each of these objects one must consider all possible distributions of 
the factors D already mentioned, which is a matter of combinatorics. In the procedure one 
must take into account the location of each D within any of the N slots available in each 
additive term and, additionally, the combinatorics associated with the choice of angles, i, i', 
for each D. 

All in all, using the quantities (2.11) and (2.10), and defining the new ones, 

r\ 1 ^ V ('S-1)! 

2"* ^ V ^ / (s - 2m + 2kj] '(s -2k- 1)!' ’ 

d I 2 

FNiA,A') = Nj2(^J^^^) {do-l)C[N-l-l,do-l] {-A')i-Ayil-Af-^-\ 

( An alternative form for F]^{A,A') is: 


N-l N-l-l 


Fr,(AA') ^ 


m 


dA 


1=0 k=0 


2^-^-i k\l\{N -k-l- 1)! {do-2k-I- 1)! 


do! 


{do + 2k-2{-l + N 
we end up with the result 


(_A'(r))(-/l(r))'(l - AM) 


7V-Z-1 


^[Cn]o = A) + r^^°-^^PM{A) , 
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but it turns out that 


^Pn{A) = {2do -2N + l)FNiA, A '), 
dr 


and therefore 


— ^ 2do-2N+lp, 

4^ ^ 2do-2N + 1 dr\ 


N 


which proves our claim. 


Appendix B: Proof of (2.7) 

Here we prove (2.7). To this effect we will consider the Noether identities for gauge theories 
applied to the Lovelock Lagrangians. These Lagrangians are generally covariant -and hence 
gauge theories- under the diffeomorphism group (coordinate reparametrizations). 

Consider a field theory with dynamics derived from a variational principle with La- 
grangian L and fields (here (j)^ means a generic field or field component). The Noether 
symmetries (satisfying that 5C is a divergence) are implemented by infinitesimal transforma¬ 
tions 

5(t>^ = , ( 6 . 2 ) 

with e“ being the infinitesimal parameters of the symmetries, and with the index a rnnning 
over the number of independent symmetries. The coefficients R^, Ra^ are functions of the 
fields and their derivatives. The gauge character of the symmetries relies in the fact that the 
functions e“ are arbitrary functions of the coordinates -and even of the fields. If \C\a are the 
functional (Enler-Lagrange) derivatives with respect to the Noether identities are written 
as 

= (6.3) 

Applying these identities to the Einstein-Hilbert Lagrangian one obtains the well known 
doubly contracted Bianchi identities for the Riemann tensor. 


For convenience, we will move momentarily to the formalism that uses the vielbein and 
the Lorentz connection as the independent fields. The vielbein extracted from the metric 
(2.1) is (we use the notation I = ^ for the flat indices) 


eg- = VA(r), 


1 

TWy 




(6.4) 


where fi and /j/ are trigonometric functions of the angles of the spheres whose specific form 
is irrelevant to ns. notice that the vielbein is “diagonal” in its indices. 

The Lovelock Lagrangian (with the cosmological constant term) in this formalism will be 
denoted Cn and its EOM are the fnnctional derivatives with respect to the vielbein and the 
Lorentz connection. The functional derivatives with respect to the connection are proportional 
to the torsion tensor and we require this tensor to vanish in order to recover onr original setting 
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where the connection was just the Levi Civita connection. Thus we need not worry about 
the EOM coming from variations of the connection. 

There are two types of gauge symmetries in the present formalism, namely diffeomor- 
phism invariance and local Lorentz invariance. The Noether identities for the local Lorentz 
invariance are empty in our case. So we will focus on diffeomorphism invariance. The vielbein 
vectors transform, under a difeomorphism 6 x^ = —e^, with the Lie derivative 


6e 


i = 


+ epdpeP 


and thus we infer from (6.2) the identifications 






^p ■ 


For the N-Lovelock Lagrangian the Noether identities (6.3) become 

- d.{[CN]iei) 44^"") = 0. (6.5) 

As an identity, (6.5) holds for any field configuration, and in particular for (6.4) (which 
implies (2.1)). When we plug (6.4) into (6.5) we can ease the notation because, the vielbein 
being “diagonal”, each contributing field is identified by a /r index: —>■ /r, we get 

4n]p R4 - dA[^N]p = 0 . ( 6 . 6 ) 

It is easy to see that the functional derivatives [Cj\f]p are proportional to those in (2.4). To 
be specific, the result is 

[Cn]p =—fi[CN]p (6.7) 

e-p 

(clearly no sum for /r). 

The only nontrivial case for (6.6) is when p = 1 (the radial variable). We obtain, using 
(2.5), (2.6) and (6.7), 


0 = [CN\pdie^ - dp{[CN\pe^) = [Tat]^( 9ie^ - 5i([£Ar]i ep) (6.8) 

= [^N]i diel + [C]si]i' 9ie| + [TatJo die^ + [Cn]i Siey — 0i(|£Ar]i ep) (6.9) 
= [£-N]i fi + [£-N]i' fi' + [-CAf]o die^ — cli(|£Ar]i) Cf (6.10) 

= —( 6 . 11 ) 

r ^/Mr) aMH 

= ^[£jv]i-5i([£^]o). (6.12) 


which is the result (2.7). 


- 18 - 








References 


[1] N. Dadhich, S. G. Ghosh and S. Jhingan, “The Lovelock gravity in the critical spacetime 
dimension,” Phys. Lett. B 711 (2012) 196 [arXiv:1202.4575 [gr-qc]]. 

[2] N. Dadhich, “Characterization of the Lovelock gravity by Bianchi derivative,” Pramana 74 
(2010) 875 [arXiv:0802.3034 [gr-qc]]. 

[3] D. G. Boulware and S. Deser, “String Generated Gravity Models,” Phys. Rev. Lett. 55 (1985) 
2656. 

[4] J. T. Wheeler, “Symmetric Solutions to the Maximally Gauss-Bonnet Extended Einstein 
Equations,” Nucl. Phys. B 273 (1986) 732. 

[5] J. T. Wheeler, “Symmetric Solutions to the Gauss-Bonnet Extended Einstein Equations,” Nucl. 
Phys. B 268 (1986) 737. 

[6] B. Whitt, “Spherically Symmetric Solutions of General Second Order Gravity,” Phys. Rev. D 
38 (1988) 3000. 

[7] M. Banados, C. Teitelboim and J. Zanelli, “Dimensionally continued black holes,” Phys. Rev. D 
49 (1994) 975 [gr-qc/9307033]. 

[8] J. Crisostomo, R. Troncoso and J. Zanelli, “Black hole scan,” Phys. Rev. D 62 (2000) 084013 
[hep-th/0003271]. 

[9] N. Dadhich, J. M. Pons and K. Prabhu, “Thermodynamical universality of the Lovelock black 
holes,” Gen. Rel. Grav. 44 (2012) 2595 [arXiv:1110.0673 [gr-qc]]. 

[10] X. O. Camanho and J. D. Edelstein, “A Lovelock black hole bestiary,” Class. Quant. Grav. 30, 
035009 (2013) [arXiv: 1103.3669 [hep-th]]. 

[11] H. Maeda, S. Willison and S. Ray, “Lovelock black holes with maximally symmetric horizons,” 
Class. Quant. Grav. 28 (2011) 165005 [arXiv: 1103.4184 [gr-qc]]. 

[12] H. Nariai, Sci. Rep. Tohoku Univ. 34 (1950) 160; 35 (1951) 62. 

[13] G. Dotti and R. J. Gleiser, “Obstructions on the horizon geometry from string theory 
corrections to Einstein gravity,” Phys. Lett. B 627 (2005) 174 [hep-th/0508118]. 

[14] C. Bogdanos, C. Charmousis, B. Gouteraux and R. Zegers, “Einstein-Gauss-Bonnet metrics: 
Black holes, black strings and a staticity theorem,” JHEP 0910 (2009) 037 [arXiv:0906.4953 
[hep-th]]. 

[15] G. Dotti, J. Oliva and R. Troncoso, “Vacuum solutions with nontrivial boundaries for the 
Einstein-Gauss-Bonnet theory,” Int. J. Mod. Phys. A 24 (2009) 1690 [arXiv:0809.4378 [hep-th]]. 

[16] G. Dotti, J. Oliva and R. Troncoso, “Static solutions with nontrivial boundaries for the 
Einstein-Gauss-Bonnet theory in vacuum,” Phys. Rev. D 82 (2010) 024002 [arXiv:1004.5287 
[hep-th]]. 

[17] H. Maeda, “Gauss-Bonnet black holes with non-constant curvature horizons,” Phys. Rev. D 81 
(2010) 124007 [arXiv:1004.0917 [gr-qc]]. 

[18] J. M. Pons and N. Dadhich, “On static black holes solutions in Einstein and 
Einstein-Gauss-Bonnet gravity with topology SO(n) x SO(n),” arXiv: 1408.6754 [gr-qc]. 


- 19 - 



[19] N. Dadhich, J. M. Pons and K. Prabhu, ”On the static Lovelock black holes,” Gen. Relativ. 
Grav. 45 (2013) 1131 [arXiv: 1201.4994 [gr-qc]]. 

[20] B. Kol, “Topology change in general relativity, and the black hole black string transition,” 
JHEP 0510, 049 (2005) [hep-th/0206220]. 

[21] N. Farhangkhah and M. H. Dehghani, “Lovelock black holes with nonmaximally symmetric 
horizons,” Phys. Rev. D 90, no. 4, 044014 (2014) [arXiv:1409.1410 [gr-qc]]. 


- 20 - 



